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We reanalyze the behavior of Friedmann-Lemaitre-Robertson-Walker cosmologies in the recently 
proposed quasidilaton massive-gravity model, and discover that the background dynamics present 
hitherto unreported features that require unexpected fine-tuning of the additional fundamental 
parameters of the theory for an observationally consistent background cosmology. We also identify 
new allowed regions in the parameters space and exclude some of the previously considered ones. 

The evolution of the mass of gravitational waves reveals non-trivial behavior, exhibiting a mass- 
squared that may be negative in the past, and that presently, while positive, is larger than the 
square of the Hubble parameter, Hq. These properties of the gravity-wave mass have the potential 
to lead to observational tests of the theory. While quasidilaton massive gravity is known to have 
issues with stability at short distances, the current analysis is a first step toward the investigation 
of the more stable extended quasidilaton massive gravity theory, with some expectation that both 
the fine-tuning of parameters and the interesting behavior of the gravity-wave mass will persist. 
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I. INTRODUCTION 

The standard cosmological model, ACDM, describes 
the acceleration of the universe by properly adjusting 
the cosmological constant A. While this simple model is 
consistent with current observational data, other models 
provide alternative explanations of this acceleration. For 
example, some models attribute the acceleration to the 
presence of a dynamical component known as dark energy 
and others to a modihcation of the gravitational 
laws on cosmological distances [IHH]- The questions will 
be to what extent it is possible to discriminate among 
the different models from observations, and whether any 
of the models are better at fitting the data than what is 
currently the most parsimonious explanation, A. 

The next generation of experiments (such as EUCLID 
[S]orDESI[in]) will provide an unprecedented amount of 
observational data. However, there is now a wide range 
of candidate theories. For instance, different modifica¬ 
tions of general relativity primarily in the infrared have 
been considered by many authors (see m for a recent re¬ 
view), and probably still others have yet to be proposed. 
Ultimately, the predictions of each candidate model must 
be confronted with data. This includes not just cosmo¬ 
logical data but data on all scales where the models make 
calculable predictions that can be tested observationally 
or experimentally. 

Within one interesting class of theories, the current 
acceleration era is associated to the presence of a mass 
term for the graviton (for a historical overview, motiva¬ 
tions and an updated description of different proposed 
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massive gravity theories, see [mn]). Here we consider 
a particular modification of general relativity known as 
quasidilaton massive gravity (QDMG), which we sum¬ 
marize in the Section |ll] This theory was proposed in 
as an extension of the dRGT theory of massive 
gravity [HI [H], and contains an additional scalar de¬ 
gree of freedom: the quasidilaton. The main motivation 
for such extension is the absence of isotropic and ho¬ 
mogeneous cosmological background solutions in dRGT 
m Indeed, it has been shown that QDMG has solu¬ 
tions with spatially flat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) background metrics [H]. Moreover, it 
has been found that (even in the absence of a cosmolog¬ 
ical constant) there are solutions for which at late times 
the metric approaches to a de Sitter metric, providing 
a plausible (self-accelerating) explanation of the acceler¬ 
ated expansion of the universe [H na [H]. The qua¬ 
sidilaton theory has three parameters more than ACDM. 

In this paper we will perform a careful analysis of the 
background cosmological evolution, taking into account 
the main goal of describing the observed expansion his¬ 
tory of the universe. While other authors have made pre¬ 
liminary investigations [Hlllg] of the background evolu¬ 
tion in QDMG, a more detailed reexamination reveals im¬ 
portant new insights. The allowed set of parameters split 
into two disconnected regions characterized by “low” and 
“high” values of a dimensionless parameter of the theory, 
w (which multiplies the kinetic term of the quasidila¬ 
ton). In the region with low values of w, while viable 
background solutions exist for a wide range of values of 
the Lagrangian parameter nominally called the graviton 
mass Trig, with rrig ^ 0 {Hq), a careful hne-tuning of the 
dimensionless constants and 0:4 is required. The per¬ 
mitted values of , 04 and nrig thus describe a very thin 
2-dim surface 771^(0:3, a4) in the {0:3,0:4,7713} parameter 
space. In the other region, the parameter rrig is con- 
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strained to be much smaller than i7o, and the larger it 
is, the narrower the 2 -dim surface of allowed and 04 . 

The paper is organized as follows. After summarizing 
the theory QDMG in Section|TTj in Section [ITII we present 
the dynamical equations. In Section |IV| we analyze the 
existence of viable de Sitter fixed-point attractors. By ex¬ 
ploring the 4-dimensional parameter space of the theory, 
in Section |V] we assess the viability of a self-accelerating 
explanation of the current expansion of the universe. An 
important outcome of our analysis is that, in order to 
reproduce an expansion history consistent with data, the 
graviton mass parameter must also be fine-tuned to a 
value that depends on other parameters of the model. 

In Section |VI| we study the evolution of the mass of 
gravitational waves for the allowed set of param¬ 

eters. We find the current value of Mqw to be gener- 
ically larger than the current Hubble constant Hq even 
when we set the graviton mass parameter rUg Hq. In 
the past (for example at redshifts relevant for the Cos¬ 
mic Microwave Background) Mqw can be either real or 
imaginary. For a conservative choice of 0 < rUg < Hq, 
\Mcw{t)\ < H{t) in the past, with IMcvyj < at 

last scattering. While this precludes the development of 
a catastrophic instability when Mqw is imaginary, never¬ 
theless potentially there could be observable cosmological 
signatures. These merit further investigation pOl - E^ . 


II. THEORY OF QUASIDILATON MASSIVE 
GRAVITY 


We consider the action for the quasidilaton theory m- 

S = Seh + 




R- 


(1) 


^P1 


+ 2mg(/l2 + <^ 3^3 + 04/34) 


where Mpi is the Planck mass and, in addition to the 
Einstein-Hilbert action Seh, a contribution charac¬ 
terizes the quasidilaton scalar field a. In addition to the 
quasidilaton kinetic term, includes three interaction 


terms: Here 

( 2 ) 

C,= l{[ICf-3[JC][IC^]+2[IC^]), (3) 

o 

£4^^([/C]4-6[/C]2[/C2] + 3[/CT + 

+ 8[/C][/C3]-6[/C4]), (4) 

with square brackets denoting a trace, and 

. (5) 

The non-dynanmical “fiducial metric” is built from four 
Stiickelberg fields 4>°‘ (a = 0, • • • ,3), 

Z/ii/ = (6) 


In the space of Stiickelberg fields, the theory enjoys the 
Poincare symmetry |I4j 

^ + c“, ^ a;)^^ , (7) 

and in addition, there is a global symmetry given by 

a^a + ao, ^ ^ 

with tTo an arbitrary constant. 

The addition of So- to the action, introduces four new 
parameters: the dimensionless kinetic coupling w, the 
graviton mass parameter rrig, and the coupling constants 
a 3 and 04 . As shown below, the cosmological solution 
depends sensitively on the values of these parameters. 


III. THE BACKGROUND COSMOLOGICAL 
EQUATIONS 


We consider a spatially flat Friedmann-Lemaitre- 


Robertson-Walker ansatz, for which 

= —N{tydt^ + a{t}'^dijdx’‘dx^ , (9) 

0° = <(’°(f), (10) 

<P^ = x\ (II) 

a = d{t). (12) 

The fiducial metric reduces to 

foo = -n[tf , fij = Sij , (13) 

where 

(14) 


he minisuperspace action for the background metric 
ad fields can now be written as 


S/V = [ dt\-3 


a3 
N 


r> w 
- 

Mpi 2N 


-\- Na ^ Tn ^{£ l 2 T 0:3/33 -I- 04/34)} 


(15) 


£2 =3(A-l)(-2 + A(l + r)), (16) 

/:3 = -(A-l)2(-4 + A(l + 3r)), (17) 

£4 = (A-I)3(-l + rA), (18) 

d we have defined 

X = , (19) 

r^^a. (20) 

Varying the action with respect to <(>°(t) leads to 


a 


G2(A) 


= a [a^G2(X)] =0, (21) 
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where G^iX) = X(1 — X)J{X), with 

J(X) = 3 + 3(l-X)a3 + (l-X)2a4. (22) 

We use time reparameterization freedom to set TV = 1. 

In summary, the independent background equations 
are: 


the constraint equation (21), or its integral 


G 2 {X) 


c 


(23) 


and 


Ha 


pm 

3M|^’ 


Pr 

3M|^’ 


"-g 

iV2 


Gi(X), 


UJ 



2 


(29) 

(30) 

(31) 

(32) 


• the Friedman equation, 
2 

2 VMpiy 

where 

. X-1 r 


3H^ = ‘^(^y + 3mlG, (X) + ^ 1 ^, 


and radiation; 


tained from Sr, 


& = MpiH 1 - 


4G2(X) \ 
XG'^iX)) ■ 


(33) 


(24) Eq. (23) implies that as a ^ oo, X ^ constant. There¬ 


fore the set of variables {ilm, ^r, X, IIa} will approach 
constants in the asymptotic future. We study the dy¬ 
namical stability of the system by means of the following 
equations: 


(34) 

(35) 

(36) 


(25) 


-2ilr ( 

'2 + 4 

matter 

dN 

r 1 



dQ\ 

-211 a 


sor ob- 

dN 

^ GiG' 


dX 

-4^ 


0 , (26) 

dN 

G '2 



H 

IP 


where a prime means derivative with respect to X ; 


• the conservation of the stress-energy tensors of 
matter, pm = —iHpm and of radiation p^ = AHpr- 


(Note that using the constraint equation (21) one can 
show that the equation obtained by taking the variation 
of Sr, with respect to a is not an independent equation.) 


IV. DE SITTER FIXED POINT ANALYSIS 


We start by investigating the future background evo¬ 
lution of the quasidilaton massive gravity model. The 
ACDM concordance model predicts the universe will ap¬ 
proach a de Sitter phase in the future. Though we do not 
know the future of the universe, we require our model to 
reproduce this prediction, consistent with recent practice 


[MUSI- 

We rewrite eq. (24) in terms of the relative energy den¬ 
sities 


1 — ^DE + (27) 

where 

rioE = , (28) 


^ can be obtained by differentiating the first Friedmann 


equation (24) to obtain 


H 

IP 


911 ™-b 1211^-b 12^ 
^2 

liio. -Lidi An —Alh_)2 
Gi “A-I-g dj l,r 

UJ 

1 - 

4 ^2 

2 

-6 


(37) 


As noted above, we focus on de Sitter fixed points, and 
require that these critical points are attractors. The de 


Sitter critical points relative to the system (34), (35) and 
(36) are given in Table |lj where 


A± = l 


3 03 ^ /9a| 

2 0:4 y 4 a| 


3 

0:4 


(38) 


To assess the stability we compute the matrix form 
of the perturbation equations linearized around each of 
the fixed points. Then, the linear asymptotic stability of 
each fixed point can be studied by analyzing the signs of 
the eigenvalues of that matrix. If the sign of the real part 
of every eigenvalues is negative, then the critical point is 
an attractor. The results are shown in the last column 
of Table HI 

This analysis indicates there are three possible late¬ 
time de Sitter fixed points. A, B and G. For each, con¬ 
straints on the parameters and 04 are obtained by 
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F. P. 

D, 

Da 

X 


Da 

Existence 

Stability 

Eigenvalues 

A 

0 

1 — uj/6 

X+ 

0 

ta/G 

0 <C 6 (j <C 6 , 0 ^ -^+; ^ ^ 

Attractor 

-4,-4,-3 

B 

0 

1 -ta/e 

X- 

0 

ca/G 

0 < ca < G, 0 < A, A e 77 

Attractor 

-4,-4,-3 

C 

0 

1 - 3uj/2 

0 

0 

3ta/2 

0 < ca < 2/3 

Attractor 

-4,-4,-3 


TABLE I. De Sitter fixed points. 



X 


FIG. 1. The blue line shows the G 2 {X) function, while the 
red points are the G 2 {Xo) = C values for {rUg = OAHq, uj = 
0.01, 03 = 0.75, 04 = 0.345}. 


requiring that > 0 and X > 0. For the point A we 
obtain 03 > 0 and 0 < 04 < 2q;§/ 3; for B, 03 < —3 and 
—3 — 303 < Qf 4 < 2q!§/ 3; for C, 04 < —6 — 403 . Notic¬ 
ing that one must insist that X > 0, we found different 
constraints than nann]. 

Consider more closely the fixed point B. Given the 03 
and 04 constraints for B, we obtain 0 < X- < X+ < 1. 
The constraint equation (231 implies that in the asymp¬ 
totic future 6*2 (X) = 0. Moreover (23 1 requires that G 2 
should be unbounded either above or below in order to 
have a past history. G 2 is a polynomial in X, so this is 
impossible if 0 < < 1 [23] as it is. Thus B cannot 

be a well-defined fixed point. Recalling that X > 0, a 
similar argument can be applied to the point C. 

The allowed {03,04} parameter region for A entails 
that 1 < X- < X+. Therefore A is the only well-defined 
de Sitter fixed point for the QDMG theory. We empha¬ 
size that our findings now differ from those of HKIS], in 
that we exclude the points B, C. 


V. COSMOLOGICAL EVOLUTION AND 
PARAMETER FIXING 

The aim of this section is to study the evolution of the 
relevant background quantities in agreement with the re¬ 
sults of the previous section and with the observed cos¬ 
mological history. That depends on the initial conditions. 


on the expansion history and on the fixed point A. By 
means of this analysis we constrain the four parameters 
of quasidilaton massive gravity: {uig^ui, 03 , 0 : 4 }. 

Given that we are dealing with the background en¬ 
ergy density evolution, we can consider neutrinos to be 
relativistic, since the value of is not negligible only in 
the radiation era when neutrinos were indeed relativistic. 
The spectrum of the CMB today is precisely measured, so 
we accurately determine U.y_o- For relativistic neutrinos, 
is proportional to Therefore we assume is 

known and we fix it by ^r,o = ^ 7.0 + = 0.0000851. 

To be definite we also fix flDE,o = 0.72, close to the 
best fit value [24]. That corresponds in ACDM to Zeq ~ 
3300. We stress that this choice will not qualitatively 
affect our conclusions on the quasidilaton massive gravity 
background evolution. 

To fix the initial conditions we require Q qe. o = ^< 7,0 + 
Ua_o, where Uo.,o and are given by (31), (32). In 
this way we obtain a 9-th order polynomial that has no 
analytical solutions. In Fig. we plot the G 2 (A) func¬ 
tion for {mg = 0.4i7o, w = 0.01, as = 0.75, 04 = 0.345}. 
The red points represent the values of Xq and G 2 (Ao) 
where the initial conditions are satified. From eq. (23) 
we obtain G = G 2 (Ao). Finally, notice that G 2 is un¬ 
bounded as A —>■ 00; this implies that for our model the 
correct past evolution of the background is allowed only 
if Ao > A+. 

The dark energy equation-of-state parameter wde,o 
is constrained by observations. To compute w for the 
QDMG model we first define the total effective equation- 
of-state parameter 


Weff- 3^2, 


and consequently 


wde 


Wr^r 

Ua -|- Dcr 


(39) 


(40) 


We must require that —1.2 < woe.o < —0.9 in agreement 
with the current limits [25] . 

The quasidilaton massive gravity model shows a partic¬ 
ular feature - ^I^e scales as matter at early times [18] . 
Indeed, from the analysis above we have 04 > 0 and 
G > 0. At early times G 2 ^ 04 A"^ = C/a'^. Therefore 
we find [26] 


noEH"^ ^ ^ m\ 


-) 

CtiJ 


3/4 


a-3 . (41) 
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FIG. 2 . (REGl): ml/HQ — oj constraints after marginalizing 
over as — 04. 



FIG. 3 . (REGl): as — 04 constraints for different ml/HQ 
values {black bold numbers). We marginalized over ljJ. The 
blue line corresponds to the boundary of the region 04 < 
2 q:|/ 3 , which is the existence condition we obtained from the 
fixed point analysis. 


It follows that, for redshifts z > 10, Ude would con¬ 
tribute to the effective matter energy density! Therefore 
^de should be negligible in the radiation era in order 
to have a viable expansion history. We demand that 

^DE{Zeq) < 0 . 01 . 

In order to identify the allowed ranges of the four pa¬ 
rameters of QDMG, the main computational obstacle 
is to find the solutions of the initial condition, namely 
^DE,o = 0.72. In principle floE.o = 0.72 could have from 
1 to 9 allowed solutions for each value of the QDMG pa¬ 
rameters. However, after enforcing all the observational 
conditions, we find that there is never more than 1 viable 
solution. We identify two disconnected allowed regions 



0.002 0.004 0.006 0.008 0.010 
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FIG. 4 . (REG 2 ): ml/Hs — uj constraints after marginalizing 
over as — 0:4. 


in the 4-dimensional space of parameters, one shows just 
low-w values (hereafter REGl) and the other one high-w 
values (hereafter REG2). 

In Fig. and we present the constraints for the 
(REGl) parameter space. After marginalizing over — 
Ui we find that uj is constrained to 0 < w < 1.2 as we 
report in Fig. On the other hand, marginalizing over 
UJ, the contour plot reported in Fig. shows that the 
as — 04 values are tightly related to the value 

{black bold numbers). Once we know two of the three 
{mg, as, 04 } parameters, the other one is determined to 
a good approximation. In other words the quasidilaton 
massive gravity theory presents a fine-tuning of the pa¬ 
rameters. 

Repeating the procedure for (REG2) we find a different 
behavior. In Fig. 0 we see that m'^/HQ < 0.008 while 
the UJ allowed interval depends on mg and it becomes 
larger as mg decreases. The 03 — 04 region is again mg 
dependent, however the dependence now is different than 
for (REGl) as we show in Fig.For (REG2) ifrUg/Tdg ^ 
0.1 then UJ 6 and 04 becomes effectively a function of 
as, so we find a fine-tuning of 2 parameters. On the other 
hand if the graviton mass is small, i.e. m'^/Hg < 0.001, 
the other parameters are no longer strongly constrained. 

Notice that we find different results than [111119]. In 
particular they allowed uj to be negative and they obtain 
UJ < 0.3. 

In our analysis we did not compute the whole expan¬ 
sion history for each point in the four dimensional pa¬ 
rameter space for practical computational reasons. As 
an illustrative example, we choose two set of allowed pa¬ 
rameters for (REGl) and (REG2) and we plot the evo¬ 
lution of the energy densities in Fig. As expected, the 
two panels, are consistent with the observed expansion 
history. 

The parameter fine-tuning we found practically re¬ 
duces from four to three the effective parameters of the 
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FIG. 5. (REG2): as — 04 constraints for different ra^/HQ values {black bold numbers). We marginalized over uj. The blue 
line corresponds to the boundary of the region 04 < 2a|/3, which is the existence condition we obtained from the fixed point 
analysis. 


quasidilaton massive gravity theory. We expect that 
studying the perturbations will further constrain the the¬ 
ory. Some of those perturbations will be unstable. 


VI. GRAVITY WAVES 

In this section we focus on the evolution of the mass 
of the gravitational waves. We consider tensor perturba¬ 
tions around the background metric solutions, 

6g^j = J ( 2^)372 exp(7fc • x) -f cc, (42) 

with = 0, and = 0. After a straightforward 

calculation, one gets the quadratic Lagrangian for 

ij,k 


Cgw = 


ITT 1 2 ' 


\h. 




TT 12 


GW I 


, (43) 


where the mass of the gravitational waves Mgw is given 

by [17] 

Mqw ~ (3 4” 30^3 -f 04 — (It 2^3 -f a4)(l -t- t)X 

+ (as -I- a4)rA^) . (44) 

We start by computing the ratio at redshift 

z = 1100, relevant for CMB, for the two disconnected 
regions (REGl) and (REG2) defined in the previous sec¬ 
tion. An exploration of the values computed reveals there 
is a minimum and maximum for each of the 

regions. The results are presented in Table jllj We see 
that a real mass as large as Mgw ~ 1Q~^H can be ob¬ 
tained, even for our conservative choice rUg < Hq. For 
both parameter regions (REGl) and (REG2), we note 
the mass can be imaginary. However, the maximum ab¬ 
solute values turn out to be much smaller than the Hub¬ 
ble rate, preventing the development of a full instability. 
It is worth noting that so far signatures in the Gosmic 
Microwave Background (CMB) due to a non-vanishing 
Mgw have been studied assuming this mass is always 
real Our results suggest that one should explore 
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Logio(1+z) 



FIG. 6. Upper panel-, expansion history for {m^/Ho = 
0.05, uj = 0.9484, as = 7.059, Q 4 = 10.63}. Lower panel: 
expansion history for {m^/Hs = 0.0005, uj = 5.4211, qs = 
10.797, 04 = 17.680} 


also the possibility of having cosmological gravitational 
waves with a small but imaginary mass at the relevant 
redshifts for CMB. 

In Fig. we plotted the evolution of the ratio 
Mqw/H^ for the parameters given in Table [n| We no¬ 
tice that at low redshifts (and in particular at z = 0) the 
mass becomes positive, and is larger than Hq, despite 
the fact that rUg < Hq. In order to assess the generality 
of this result we computed the ratio Mq^/Hq at z = 0 
varying the parameters in the two disconnected allowed 
regions of the 4-dimensional space, and we obtained its 
maximum and minimum value. The results are shown in 
Table [Till We see that Mqw is larger than Hq, even for 
values of rrig ~ I0“^i7o, and it is up to a factor of ~ 5 
larger than Hq for rrig < Hq. The existence of a mini¬ 
mum value of Mqw that is larger than Hq is remarkable, 
since this represents a motivated observational thresh¬ 
old. That is, if one could constrain Mqw to be smaller 
than Hq one would be able to rule out a self-accelerating 
explanation of the current acceleration of the universe 
within the QDMG theory. It would be interesting to see 
whether an analogous result holds for other theories that 


also aim to provide a self-accelerating explanation. Un¬ 
fortunately, current experiments are still far from probing 
Mow ~ Hq [^. Moreover, the upper limits one can ob¬ 
tain are in general model dependent, since they are based 
on assumptions involving different scales of the theory. 
This represents a challenge for both theory and observa¬ 
tions, and highlights the ongoing importance of working 
out predictions within the framework of specific models 
of modified gravity. 


0 

0 

t-H 

II 

(REGl) 

Max. Min. 

(REG2) 

Max. Min. 


2.5 X 10“'' 

-9.2 X 10“® 

1.5 X 10"® 

-8.5 X 10"® 

ml/Hi 

1 

1 

lO"'^ 

8.45 X 10“® 

OJ 

1.03 

0.32 

5.99 

5.95 

as 

6.60 

6.80 

5.83 

3.56 

a 4 

19.61 

19.87 

20 

8.44 

Xo 

1.82 

1.85 

1.73 

1.94 


TABLE II. Maximum and minimum values of at 

2 = 1100 and the corresponding values of the parameters 
{rrig, a;, 03 , 04 }. For completeness, the value of X at 2 = 0 is 
also presented. 



FIG. 7. Evolution of the mass-squared of the gravitational 
waves in units of the Hubble rate for the set of parameters 
given in Table|^ Solid (Dashed) lines correspond to the values 
of parameters for which we found the maximum (minimum) 
value of Mq^y/H^ at 2 = 1100. 


VII. CONCLUSIONS 

The combination of General Relativity and the Stan¬ 
dard Model of particle physics are demonstrably and re¬ 
markably successful descriptions of the world on scales 
up to and including the solar system. On larger scales, 
there is a need either to modify the theory of gravity 



































z = 0 

(REGl) 
Max. Min. 

(REG2) 
Max. Min. 

Mhw/Hi 

25 

2.5 

18 

5 

ml! Hi 

1 

10 “'‘ 

lO""^ 

10 “'‘ 

UJ 

0.08 

0.08 

6 

5.05 

Q3 

0.86 

3.46 

5.83 

1.70 

CX4 

0.48 

0.55 

20 

0.49 

Xo 

4.67 

19.59 

1.73 

10.99 


TABLE III. Maximum and minimum values of 

at 2 = 0 , the corresponding values of the parameters 

{mg,uj, as, 04 }, and the value of X at z = 0 . 

or to introduce new forms of dark matter and dark en¬ 
ergy. The most parsimonious solution would be to iden¬ 
tify candidates for the latter in the Standard Model, and 
such candidates may exist for dark matter (see for exam¬ 
ple [29l [30] ) and evade existing constraints m, although 
the phenomenological successes of MOND (see for exam¬ 
ple [32]) cannot be entirely dismissed as an indication of 
the need to modify gravity on galactic scales. For the 
observed cosmic acceleration, the situation is even less 
clear. A cosmological constant is the canonical explana¬ 
tion, but despite decades of attempts has as yet no clear 
explanation in the Standard Model. The need for ob¬ 
servational probes of possible dark energy and modified 
gravity explanations is thus paramount. 

One possibility would be to develop some general 
phenomenological classification of possible deviations of 
gravity from GR. The Parametrized post-Newtonian ap¬ 
proach is one such program, in the context of almost- 
Schwarzschild backgrounds. Such generic approaches 
have also been attempted in the cosmological context 
(eg. [33|[34]). However, in the context of a highly non¬ 
linear theory such as GR, the observational consequences 
of small theoretical departures from GR can be quite 
ideosyncratic. While phenomenological parametrization 
of observables may be convenient, and even useful, they 
may not capture (or may capture poorly) the specific 
phenomena or behaviour that result from actual mod¬ 
els. Careful examination of specific individual models 
can therefore be both instructive and essential. 

In this paper, we have studied the (homogeneous) cos¬ 
mological solutions of quasidilaton massive gravity. A 
study of the linear perturbations around the asymp¬ 
totic self-accelerated cosmological solution of this theory 
(which corresponds to a De Sitter background metric) 
has been done in jSS] [35] . These studies have revealed 
that the kinetic term of one of the scalar perturbations 
becomes negative for short wavelengths, indicating that 
the theory may have a ghost instability that shows up at 
short distances. This is indeed the case at linear level. 
Several authors jST] [38] have therefore extended the the¬ 
ory by allowing for a new coupling, which can be prop¬ 
erly adjusted to make the scalar sector stable at linear 
level. This extended quasidilaton massive gravity theory 


(EQDMG), has been considered by other authors [531441] . 

Although this current reconsideration of the back¬ 
ground cosmological solutions of QDMG was performed 
as a first step for a full analysis of the EQDMG, it re¬ 
vealed important attributes of the QDMG cosmology, 
which we expect to carry over qualitatively or in detail to 
EQDMG. The first is that observationally viable QDMG 
cosmologies require fine-tuning of parameters. In par¬ 
ticular, the allowed values of the graviton mass parame¬ 
ter, irig is a tightly constrained function of the coupling 
constants and 04 , with only a very narrow tolerance 
around a central value 7713 ( 03 , 04 ). This fine-tuning, and 
the precise value of 7773(03,04), is dictated by observa¬ 
tional constraints on the dark energy properties. 

The second observation is that some small (but pos¬ 
sibly non-negligible) fraction of what manifests as floE 
(i.e. p/p ~ — 1 ) today, was (i.e. p/p ~ 0 ) in the 
past. The transition from one equation of state to the 
other was sudden and probably not well-captured by a 
linear parametrization of w(z). The expected difference 
between at high redshift (as measured in the GMB) 
and at low redshift (as measured, say in large scale struc¬ 
ture), could be the source of recently noted tensions in 
different determinations of fim [HI US]. While the de¬ 
tails of these behaviors of the background cosmomology 
are likely to be altered in EQDMG, it is plausible that 
these qualitative features are robust. 

We have also analyzed the phenomenology of the gravi¬ 
ton. The governing equations for the graviton mass Mow 
(which is not equal to the graviton mass parameter rrig) 
are the same in QDMG and EQDMG. We therefore ex¬ 
pect to gain useful insights for the extended model pro¬ 
vided that the background solutions do not depend sen¬ 
sitively on the new parameter ao- of the extended theory. 
We find that the graviton mass-squared typically is neg¬ 
ative at redshifts well above z = 1, indicating an insta¬ 
bility. This includes redshifts z ~ 10^ — 10^ where such 
physics may well imprint itself on the GMB. At any given 
time «C , so we do not expect the instability 

to lead to many e-foldings of growth. Nevertheless, if 
this persists in EQDMG, it may be another opportunity 
to see evidence of modified gravity in GMB observations. 

Regarding vector perturbations, according to equa¬ 
tions (4.16) and (4.17) of [5T] . the square of the speed 
of propagation, Cy, can be recast as Cy = 
where the absence of ghost instability is guaranteed pro¬ 
vided Ky > 0 [33]. Notice in particular that when 
becomes negative, the absence of ghost instability implies 
Cy becomes also negative. Therefore, we expect that a 
detail analysis of the perturbations will further reduce 
the region of allowed parameters. 

In a future work, we will therefore extend our analy¬ 
sis to the EQDMG theory, taking into account the con¬ 
strains from the study of the perturbations, anticipating 
hopefully that these observable effects will indeed persist. 
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